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Abstract 

We investigate consequences of space non-commutativity in quan- 
tum mechanics of the hydrogen atom. We introduce rotationahy in- 

^ 3 

variant noncommutative space Rg - an analog of the hydrogen atom 
(ff-atom) configuration space Rg = R'^\{0}. The space Rg is gener- 
ated by noncommutative coordinates reahzed as operators in an aux- 
ihary (Fock) space J-. We introduce the Hilbert space % of wave 
functions V formed by properly weighted Hilbert-Schmidt operators 
in F. Finally, we define an analog of the H-&iom. Hamiltonian in Rg 
and explicitly determine the bound state energies and the corre- 
sponding eigenstates i^nj-m- '^^^ Coulomb scattering problem in Rg is 
under study. 



1 Introduction 



Basic ideas of non-commutative geometry have been developed in and, in a 
form of matrix geometry, in [2j . The main apphcations have been considered 

• in the area of quantum quantum field theory in order to understand, or 
even to remove, UV singularities, and 

• eventually to formulate a proper base for the quantum gravity. 

The analysis performed in [5] led to the conclusion that quantum vac- 
uum fluctuations and Einstein gravity could create (micro)black holes which 
prevent localization of space-time points. Mathematically this requires non- 
commutative coordinates x'^ in space-time satisfying specific uncertainty re- 
lations. The simplest set of operators representing in an auxiliary 
Hilbert space should satisfy Heisenberg-Moyal commutation relations 

[x^^x"] = iO'"', /i,z/ = 0,1,2,3, (1) 

where 6'^'^ are given numerical constants that specify the non-commutativity 
of the space-time in question. 

Later in |1] it was shown that field theories in NC spaces with ([T]) can 
emerge as effective low energy limits of string theories. These results sup- 
ported a vivid development of non- commutative QFT. However, such models 
contain various unpleasant and unwanted features. The divergences are not 
removed, on the contrary, UV-IR mixing appears, [3] . The Lorentz invariance 
is broken down to SO (2) x 5*0 (1, 1), but even this is sufficient to prove the 
classical CPT and Spin-statisics theorems, [B]. This was not accidental and 
led to the twisted Poincare reinterpretation of NC space-time symmetries. 

However, it could be interesting to reverse the approach. Not to use the 
NC geometry to improve the foundation of QFT, what is a very compli- 
cated task, but to test the effect of non-commutativity of the space on the 
deformation of the well-defined quantum mechanics (QM): 

• Various QM systems have been investigated in 3D space with Heisen- 
berg-Moyal commutation relations = i6'^\ i,j = 1,2,3, e.g. harmonic 
oscillator, Aharonov-Bohm effect. Coulomb problem, see |H], [S]- However, in 
such 3D NC space the rotational symmetry is violated and there are systems, 
such as if-atom, that are tightly related to the rotational symmetry. 

• The rotational symmetry survives in 2D Heisenberg-Moyal space with 
NC coordinates Xi,X2 satisfying the J-" commutation relations [a;i,X2] = iO 
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in an auxiliary Hilbert space. In [10] a planar spherical well was described 
in detail: 

(i) First, the Hilbert space l-L of operator wave functions if) = iIj{xi,X2) 
was defined; 

(ii) Further, the Hamiltonian was defined as an operator acting inTi. It 
was nice to see how the persisted rotational symmetry helps to solve exactly 
the problem in question. 

The presented list of references is incomplete and we apologize for that. 
We restricted ourselves to those which initiated progress or are close to our 
approach. 

Our aim is to extend this scheme to the QM problems with rotationally 
symmetric potentials V{r) in the configuration space = \ {0}. We 
restrict ourselves to the Coulomb potential which, in the usual (commutative) 
setting, is a solution of the Poisson equation finite at infinity: 

AV{r) = ^ V{r) = + q^. (2) 

For if-atom q, in a Gaussian system of units, is a square of electric charge 
e^, and we put the inessential constant qo = 0. In this case we are dealing 
with Schrodinger equation 

- — AV^(x) - — ^(x) = Eip{x), r = |x| > (3) 

in the Hilbert space Ti specified by the norm 

I = J ci3x|^(x)p. (4) 

Expressing the wave function as 

V>(x) = R^{r)Hj^{x), H,U^) ~ r^Y.U^,^), (5) 



and putting a = 2me^/ ti^ and k = \/—2mE / h, we obtain the radial Schrodin- 
ger equation: 

r i?J(r) + 2(j + l)R'j{r) + aRj{r) = K^r Rj{r) . (6) 

Its solution is 

R,{r) = e-''^^ F (j + 1 - ^, 2j + 2; 2/s:r) , (7) 



where -F(a, c; x) is the confluent hypergeometric function. 

For bound states E < 0, i.e. real-valued k, the solution should have a 
finite norm in Tio- This is the case when the first argument of the degenerated 
hypergeometric function is zero or a negative integer, and this determines the 
discrete energy eigenvalues: 

n = j + j+2, ... ^ e^ = -—kI = -——. (8) 



The hypergeometric function then reduces to a polynomial of degree n—j — 1 
in the variable x = 2rK: 

F(ni-n2n2-x)=^Tc^t^ = + 



k) ■ 

(9) 

In this paper we extend the QM solution ©-(IH]) of the Coulomb problem 
in Rq to the non-commutative rotationally invariant space. In Section 2 we 
define Rq , a rotationally invariant NC generalization of the configuration 
space Rq, we introduce the generators of rotations and the NC analog of 
their eigenfunctions. In Section 3 we define Hilbert space H - the NC analog 
of T-L, and we introduce the NC analog of the Coulomb problem Hamiltonian 
acting in Ti. In Section 4 we exactly solve the corresponding NC analog of the 
Schrodinger equation. Last Section 5 contains conclusions and perspectives. 



2 The noncommutative space Rq 

In this section we define the noncommutative space Rg, possessing full ro- 
tational invariance, as a sequence of fuzzy spheres introduced, in various 
contexts, in [H]. Different fuzzy spheres are related in such a way that at 
large distances we recover space Rq with the usual flat geometry. A similar 
construction of a 3D noncommutative space, as a sequence of fuzzy spheres, 
was proposed in [12]. However, various fuzzy spheres are related to each 
other differently (not leading to the flat 3D geometry at large distances). 

We realize the noncommutative coordinates in Rq in terms of 2 pairs of 
boson annihilation and creation operators d^, aj,, a = 1,2, satisfying the 
following commutation relations, see |13] : 

[a„,ajj] = 5ai3, [da,d,3] = [a),,ajj] = 0. (10) 
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They act in an auxiliary Fock space T spanned by normalized vectors 



= ^ ^^^A^ |0). (11) 



Here, 1 0) = |0, 0) denotes the normalized vacuum state: Si |0) = 0-2 |0) = 

Cj, j = 1,2,3, in the space Rq 



The noncommutative coordinates x,-, j = 1,2,3, in the space Rq are 



given as 

= Xa'^aja = Xcripa^^ap, j = 1,2,3, (12) 

where A is a universal length parameter. The coordinates xj satisfy rotation- 
ally invariant commutation rules: 

[xi,Xj] = 2iXeijkXk, [xi,g] = 0, (13) 

where g = XN, and N = d = dl^da- 

The operator that approximates the Euclidean distance from the origin 
in an optimal way is f = A (A^ + 1) = ^ + A, and not g. Namely, it holds 
f'^—x'j = A^, whereas g'^—x'j = o(A). In Section 3 we give a strong argument 
supporting the exceptional role of f. 

Let us consider the linear space of normal ordered polynomials containing 
the same number of creation and annihilation operators: 

^ = Cm,n.,n.n, {al^ (4)"^^ i^lT' («2)"^ , (M) 

where the summation is finite over nonnegative integers satisfying mi + m2 = 
rii + n2. In this space we define generators of rotations Lj, j = 1,2,3, as 
follows 

L,^ = '-[d+a,d,% J = 1,2,3, (15) 

obeying proper commutation relations 

[Li,Lj]'ip = {LiLj - LjLi)'4) = ieijkLkip . (16) 

With respect to the rotations (|T5|) the doublet of annihilation (creation) 
operators transforms as spinor (conjugated spinor), whereas the triplet of 
NC coordinates as vector 



L,at .l.j„at. L.*,^. ...... 



The standard eigenfunctions ipjm, j = 0, 1, 2, . . . , , m = —j, . . . , +j, 
satisfying 

L'ii^jm = j{j + l)ijjm, Ls^jm = m^jm, (17) 

are given by the formula 

^ ^ „ ^ MrMM]!! <Sd^ (18, 

^ mi\m2\ ni! 71,2! 

with the summation over all nonnegative integers satisfying mi + m2 = 
ni + 71-2 = J, — 1712 — Til + n2 = 2m. Thus ijjjm = when restricted to 
the subspaces J-V = {\ni,n2) \ rii + n2 = N} with N < j. For any fixed 
: Rj{g) ■ equation (fTSj) defines a representation space for a unitary irreducible 
representation with spin j. 

The symbol : Rj{p) '■ represents a normal ordered analytic function in the 
operator g: 

: RM ■■ = Y,< = Y.^^^' (19) 

The last equality follows from the equation 

iV' 

:N^: |ni,n2) = l^i.^2), iV = m + ^2 (20) 

(which can be proved by induction in k). Since, : N'^ : \ni,n2) = for 
k > 111 + n2, the summation in f|T9|) is effectively restricted to < on any 
subspace J-at. 



3 Quantum mechanics in space Rq 

Let H denote the Hilbert space generated by functions (fTSjl with weighted 
Hilbert-Schmidt norm 

= 47rA^Tr[(iV+l)^1'^] = 47r Tr[f ^"f ^] , f = A(iV + l). (21) 

The rotationally invariant weight w{f) = AirX'^f is determined by the 
requirement that a ball in Rq with radius r = A(A^ + 1) should possess a 
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standard volume in the limit r — t- oo. The projector on the subspace 
v^o © • • • © J^N corresponds to the characteristic functions of a ball with the 
radius r = A(A^ + 1). Therefore, the volume of the ball in question in Rq is 



4vr 



Vr = 47rA^Tr[(iV+ l)Pjv] = 47rA3^(n + l)' = — r3 + o(A). (22) 

o 

n=0 

Thus, the chosen weight w{r) = Atc X'^f possesses the desired property. 

Note: The weighted trace Tr[w(f) ...] with w{f) = 4tt X"^ f at large dis- 
tances goes over to the usual volume integral J d^x ... . The 3D noncommu- 
tative space proposed in [12] corresponds to the choice w{f) = const and at 
large distances does not correspond to the flat space Rq. 

The generators of rotations Lj, j = 1,2,3, are hermitian (self-adjoint) 

operators in "H, and consequently, the two operators '^jm and ^j'm', with 
arbitrary factors : Rj{g) : and : Rj'{g) :, are in orthogonal. It is sufficient 
to calculate = ll^iilP (this equality follows from the rotational 

invariance of the norm in question): 

oo N 

W^jmf = 47rA=^ 5Z + {n,N-n\{N + l)^^^,j\n,N-n) , (23) 

N=j n=0 

We benefit from the fact that has a simple form 

% = -^(^ly --^M-- i-^2y (24) 

The matrix element we need to calculate is 

{n,N — n\ (0-2)'' • Rjio) '■ ^1 {(^{y • Rjio) • Q-2 \n,N — n) 
{n + j)\{N-n)\ 1^ 2 



T^jiN-j)\\ (25) 



n\ {N — j — n) 
where 

TZjiN) = {n,N-n\: Rj{g) : \n, N - n) (26) 

(the expression on the r.h.s. is n - independent). Inserting (p^ . ( |25|) into 
(!23|) and using the identity 

^ /^n + j\ [N -n\ fN + j + 1 

n=0 



J J \ J J V 2j + l 
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we obtain 

This expression represents, up to an eventual normalization, the square of a 
norm of the radial part of the operator wave function. 

Now we are ready to define the Coulomb problem Hamiltonian in the 
noncommutative case. 

The kinetic term in Rq. In the first we define the NC analog of the 
Laplacian in Rq as follows: 

Aa ^ = - ^ [al n = - ^2^^^ I) • ^28) 

This choice is motivated by the following facts: 

• A double commutator is an analog of a second order differential opera- 
tor; 

• The factor guarantees that the operator Aa is hermitian (self- 
adjoint) in Ti, and finally, 

• The factor A~^, or respectively, guarantees the correct physical 
dimension of and its non-trivial commutative limit. 

Calculating the action of f l3T|) on ipj^ given in f|T8l) we can check whether 
the postulate ( 13T|) is a reasonable choice. The corresponding formula is de- 
rived in Appendix A: 



\mi I ^m2 



mil 77121 

(—n ^"2 

X : [gR"{g) + 2(j + 1) R'{g)] : ' \ ' . (29) 

Uil n2i 

Here the symbols R'{g) and R"{g) are defined as: 

oo oo 

fc=0 k=l 

oo 

i^\q) = J]Mfc-i)4^'-'- (30) 

k=2 
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Thus, the prime corresponds exactly to the usual derivative dg. We see 
that the angular dependence in (12^ remains untouched, since is rotation 
invariant. In the commutative limit A — )■ formally ^ — )■ r, and we see that 
( 129|) guarantees that A a reduces just to the standard Laplacian. 

Based on that, we postulate the kinetic term of the Hamiltonian as follows 



The potential term in Rq is defined by left multiplication of the operator 
wave function by potential V = V{x): ^ i— )■ V\E'. The potential is central 
if = V{N), or equivalently, V = V{r). 



4 The Coulomb problem in 

In the commutative case the Coulomb potential is a radial solution of the 
equation (j2]) finite at infinity. Due to our choice of the noncommutative 
Laplacian Aa the equivalent equation in Rq is 

[al [a,,y(iV)]] = 0. 

It can be rewritten as a simple recurrent relation 

{N + 2)V{N + 1) ~ {N + 1)V{N)] = {N + 1)[V{N) - NV{N-1). (32) 

Putting 

(M + 1) V{M) - MV{M - 1) = go , ^(0) = ?o - I , 

and summing up the first equation over M = 1, . . . A^, we obtain the general 
solution: 

where, q and go are arbitrary constants (A is introduced for convenience). 
For if-atom g = and we put go = 0. We see that the dependence of 
the NC Coulomb potential is inevitable. 
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Thus, the noncommutative analog of the Schrodinger equation with the 
Coulomb potential in Rq is 

^^[aLK,^]]-|^ = ^^ ^ j[al[a^,^-a^ = -K^f^, (34) 



where a = 2me^/fl^ and n = \/ —2mE / h. The nocommutative corrections 
comming from f\E' are calculated in Appendix A: 

r^,^ = E ^^^r^ ■.[ig+Xj + X)R, + XgR'^]: <lzMl , (35) 

mi! 7712! ni! 712! 

where Rj = Rj{g) and similarly for derivatives, e.g., R'j = R'j{g). The 
angular dependence is again untouched as the multiplication by f represents 
a rotation invariant operator. 

Inserting f l29|) and fl35|) into fl34|) we obtain the NC analog of radial 
Schrodinger equation: 

: [gR'^ + 2{j + l)R'j + a Rj] : = : [gRj + A((j + 1)^, + gR'^)] : . (36) 

The term proportional to A represents the NC correction. For A — )■ this 
NC radial Schrodinger equation reduces to the standard radial Schrodinger 
equation. 

We associate the following ordinary differential equation to the mentioned 
operator radial Schrodinger equation (l36l) : 



gW^ + 2(j + 1) 7^; + aTZj = [gUj + A((j + 1) 7^, + gTZ'^)] . (37) 

If the function TZj = TZjig) = Tlik^kS^ solves the associated ordinary dif- 
ferential equation ( 1371) . then 

R, = : KM : = Y.^"^ : g' : = Y^c^.X' . (38) 

solves the operator radial equation Moreover, the operator function 

Rj = : TZjig) : possesses a finite norm in % provided the function TZj = 
TZj{g) has finite norm in (since the norm ( 12T|) asymptotically reduces to 
the usual QM norm). 

The solution TZj of the associated radial Schrodinger equation fl37|) is given 
similarly as in the standard Coulomb problem in ([7]), but with particularly 
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scaled q dependence in the exponent and in the argument of the confluent 
hypergeometric function: 



= e^'^' F (j + 1 - 2j + 2 ; 2QKd) . (39) 

The dimensionless quantities h and d given as (see Appendix B) 

h = y^l + T]^ - r], d = + = ^ (40) 

specify the NC corrections to the usual radial function in ([7]). They en- 
ter b and d via the parameter rj = Xk/2. So they vanish not only in the 
commutative limit A — )■ , but also for — )■ 0. 

We shall restrict ourselves to the determination of the bound states spec- 
tra with E < and k > 0. In this case TZj{g) should be normalizable. This 
is ensured if the first argument of the confluent hypergeometric function is 
zero or negative integer, what determines the discrete energy eigenvalues 
(remember that d is K-dependent, see fHU]) ): 

a 

-- n = J + 1, J + 2, ... 



2dnKji 



" = '"h^""- = ~2^ l + ^AVagn2' ^ ^ 

where ao = 2/a = h'^/me'^ = 5.29 x 10~^^m is the Bohr radius. The first 
factor in is just the standard bound state energy of the Coulomb problem, 
whereas the second one represents the noncommutative correction. The NC 
corrections in the limit A/n — )• 0, i.e., in the commutative limit A — )■ 0, or for 
fixed A, in the quasi-classical limit — )■ oo for highly excited states. 

The solutions "^nim operator equation f l5B]) corresponding to the 

energy E^ is 



(jm) 



where Nz!^ denotes the normalization factor and 



T^njig) = e-'" '^"^ F{j + l-n, 2j + 2, 2^/s:„ rf„) . (43) 
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Here, the parameters 6„ a dn are given by (gni) with K„ = ^J-2mElin. 
From the equation fl2UI) it follows directly 

■.e-^'^Q':= (1-aA)^- ttt 7-— — = (1 - aA)^/^ — . (44) 



{l-a\f{N-k)\ (1-Aa) 



This allows as to express the normal ordered form of the radial part of oper- 
ator wave function: 



: UnM := {l-hn KnXy'^ : F (^j + I - 2j + 2, 



fc=0 



with the coefficients c^^- given in dH]). Inserting 

n„,W ^ iv- (1 - '-^ ..A)" <, (t?^ . (46) 

into ( 127|1 the normalization constant N'^^ can be determined (we skip its cal- 
culation). In commutative limit the factor in front of the sum gives the usual 
exponential damping factor in ([7]). However, as the argument of the polyno- 
mial becomes scaled due to NC corrections, the separation of the asymptotic 
factor from the polynomial part is not perfect. 



5 Conclusions 

We carefully defined the NC rotationally invariant analog of the QM con- 
figuration space and the Hilbert space of operator wave functions in NC 
configuration space. The central point of our construction was the definition 
of Aa the NC analog of Laplacian, supplemented by a consequent defini- 
tion of the weighted Hilbert-Schmidt norm and a definition of the Coulomb 
potential satisfying NC Poisson equation. 

With this input this Hilbert space we introduced the NC analog of H- 
atom Hamiltonian and explicitly determined the bound-state energies 
and corresponding eigenstates ipnjm (see equations ( HTj) . fH2]) and fH5]) ). 
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We found that the discrete parameters n, j, m have the same meaning and 
range as in the standard (commutative) Coulomb problem, and moreover, the 
bound-state energies and eigenstates possess a smooth commutative limit 
A — >• 0. This paper does not deal with the case of the scattering in the NC 
configuration space - this will be discussed elsewhere. 

The noncommutativity parameter A is not fixed within our model. How- 
ever, it can be estimated by some other physical requirement. For example, 
one can postulate, as was done in early days of modern physics, that the rest 
energy mc^ of electron is equal to the electrostatic energy of its Coulomb 



field. In Rq this means: 



47rA^ 

mc^ = -—Tt[{N + 1)E^] (47) 
Svr 



where 



Ej = . X,-, (48) 

' A3 N{N+l){N + 2) ' ^ ^ 



is the NC electric field strength corresponding to NC Coulomb potential $ 
which was discussed in Section 4 (the details will be published, see |15j). 

We stress that in the NC case the electrostatic energy of electron, deter- 
mined by the trace in (HTj), is finite (no cut-off at short distance is needed). 
A straightforward calculation of the trace in P7|) gives the relation: 

mc^ = ^ A = = Ao. (49) 

8 A 8mc^ ^ ^ 

This Ao is fraction of the classical radius of electron ro = e^/mc^: Aq = 
1.06 X 10~^^m = 1.06 /m (the coincidence with the proton radius is purely 
accidental). 

The NC corrections to the H-atom energy levels given in (14T!) are of of 
order (Ao/ao) = (9/64) ~ 4 x 10"-^^ (here ao ~ 1/137 is fine structure 
constant). Such tiny corrections to energy levels are beyond any experimen- 
tal evidence. Moreover, at Ao ~ 1 fm relativistic and QFT effects become 
essential. 

Our investigation indicates that the noncommutativity of the configura- 
tion space is fully consistent with the general QM axioms, at least for the 
if-atom bound states. However, a more detailed analysis of the Coulomb 
problem in Rq would be a desirable dealing, e.g, with the following aspects: 
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• Coulomb scattering problem, dyon problem (electron in the electric 
point charge and magnetic monopole field), Pauli if-atom (non-relativistic 
spin); 

• Coulomb problem in Rq and its dynamical symmetry, QM supersym- 
metry and integrability of the Coulomb system. 

Besides non-relativistic H-atom, there are other systems that would be 
interesting to investigate within NC configuration space Ro=, e.g., Dirac H- 
atom (relativistic invariance?), non-Abelian monopoles, or spherical black- 
holes. 

Acknowledgements: The authors would like to thank to M. Chaichian, 
M. M. Sheikh- Jabbari and A. Tureanu for valuable comments. This work 
was supported by project VEGA 1/100809/1. 

Appendix A 

Here we prove two formulas, fl29l) and fl35|) . we need for the calculation of 
the NC Coulomb Hamiltonian. Below we skip indices j and m. 
a. Let us begin with the ( 129|) : 

[ai, [a„, ^jm]] = 



+ 
+ 
+ 



where R = J^'kLo '4^^' i-^-' '4 ~ '4- we shall use the following 



'"1 ( n ' V"'2 



a [a 



{jm) 



nil'. 1^2- 



4)"" (4) 

mi! 7712! 



m2 



ni\ 



ni\ 77,2! 



^ mi! m2! ni\ 712! " 

{jm) 



(a;)- (air- ,„ „ A „ <■ (-02)"" 



(jm) 



mil m2! 



[al, [a„, : ^ :]] — 



ni! 77,2! 

77i! 772! 
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commutation relations 

[al-.N^-] =-kai -N^-^: [al:R:] = -al : dj^R : , 

[a„, : N'' :] = k : iV'^"^ : a« ^ [a«, : R :] = : dj^R : a« , (51) 

where dj^ denotes the derivatives with respect to N: dj^R = J2T=i ^ ^ iV^~^ 
It is easy to see that the second hne in (!50l) vanish, and the the first and 
third hne give the same contribution 

^ m,\m^\ ^ ^ • "^^^ ' ^^^^ 

(jm) 

From fl5T]) the double commutator [aj^, [cia, '■ R '■]] follows directly, and this 
gives the value of the third line in f l50|) 

Introducing parameter A and switching to the derivatives with respect to g 
the last two equations yields ( l29l) . 

b. Proof of (1351) is straightforward. From equation ( l20l) it follows easily 

N ■.N'' : = : iV'^+i :+A;:iV'=: ^ N : R : = : N R : + : N d^R : . (54) 

This relation gives directly 

mi! ma! ni! ^2! 

= J] ... [(iV + j + 1) :^:] ... 

= J2 ... :[{N + j + l)R + Ndj^R]: ... , (55) 

where we have replaced both untouched factors containing annihilation and 
creation operators by dots. Introducing parameter A again and switching to 
the derivatives with respect to g we recover fl55|) . 



15 



Appendix B 

It is known that the solution R = R{g) of equation 

gR" + {aiQ + bi)R' + (aa^ + 62)^? = 



can be expressed in terms of a confluent hypergeometric function 
F{a, c; x) 



a X aia + 1) x"^ 



cl! c(c+ 1)2! 
The formula, given e.g. in |T3], reads 

R = e^(^-"^)^F(a,c;-D^), 
where D is determined by = a\ — 4a2 and 



D 



{D - ai) hi + 62 



hi. 



In our case 



ai 



what gives 



R[g) = ei(^+^'^^)^ F (j + 1 + ^, 2j + 2; -D g) . 



(56) 



(57) 



(58) 



-Ak^ 61 = 2j + 2, aa = -k^ 62 = a - (j + 1) A /^^ (59) 



(60) 



There are two solutions R±{g) depending on the sign oi D = ±2k a/1 + rj'^, 
T] = I A K. However, R+{g) = R-{g) due the Kummer relation 



F(a,c:x) = F(c — a,c: —x) 



We took the solution fH3l) of the associated radial NC Schrodinger equation 
fl39|) with negative D as it is more convenient for the determination of bound 
states. 
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